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ABSTRACT 


An investigation of the dynamic stability of a typical fixed wing 

VTOL aircraft operating in the transition region of flight was conducted 
by means of analog simulation. The Ryan Aeronautical Company XV-5A 
Lift Fan Research Aircraft was used as a basis for the analysis. Wind 
tunnel tests on a simple 1:24 scale model of the aircraft provided static 
aerodynamic data throughout the angle-of-attack range -909 < ox, < 
+90°. Non-dimensional stability derivative data was extended through- 
out the angle-of-attack range. Analog simulation at relative velocities 
of 100, 60, and 20 knots was conducted using dimensional parameters. 
Recordings of the simulation runs provided period, frequency, damping, 
and dutch roll parameter information which was analyzed to determine 
areas of flight instability. The variations of these parameters with both 
angle-of-attack and velocity is discussed. The need for automatic 


stabilization is indicated. 
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1.0 INTRODUCTION 


The versatility of a fixed wing V TOL aircraft will be enhanced 
if it has the capability of flying at arbitrary combinations of vertical 
and horizontal speeds in the transition from hover to conventional for- 
ward flight and vice versa. This implies flight at arbitrary angle-of- 
attack in the range - 90°S مے‎ $+90° and at arbitrary relative air- 
speeds. Throughout this flight envelope wide variations of aerodynamic 
characteristics can be anticipated. The dynamic behavior of the air- 
craft, therefore, can,be expected to change radically. 

The purpose of this analysis is to investigate the dynamic 
stability of a typical fixed wing VTOL aircraft in this transition region 
of flight throughout both the angle-of-attack range -909 Se«,, € 4909 and 
the speed range 20€ V = 100 knots. The dynamic response of such an 
aircraft after being perturbed from a reference flight condition is 
analyzed in an effort to establish the effects of the varying aerodynamic 
forces and moments on the behavior of the aircraít. Areas of insta- 
bility as well as the degree of stability of instability as functions of 
both speed and angle-of-attack are obtained. The feasibility of manual 
control in the high transition speed - high angle-of-attack region of flight 
is considered and the necessity for automatic contrcl is indicated. 

The aircraft chosen as a basis for the analysis is the Ryan 
Aeronautical Company XV-5A. This aircraft, built as a research 


aircraft for the U.S. Army, employs the General Electric X-353-5 








Life Fan Propulsion system. The aircraft is capable of both conven- 
tional flight at high subsonic speed and VTOL operations. It is not the 
intent of this analysis to define the dynamic behavior of the XV-5A per 
se but to indicate the dynamic behavior of a typical VTOL aircraft 
using the XV-5A as an example. For.,this reason, although the XV- 
5A is equipped with manual controls, only stick fixed stability was 
considered, 

The analysis required extrapolation of stability derivative data 
throughout the angle-of-attack regime. To this end, static wind tunnel 
tests were performed to obtain the basic lift, drag, and pitching mo- 
ment data. The longitudinal and lateral-directional stability derivative 
data was then extended by means of theoretical analysis, empirical 
relations, and estimates to cover the specified angle-of-attack range. 
Finally, an analog simulation study, using dimensional parameters, 
was conducted over the ranges of angle-of-attack and flight velocities 
indicated previously. 

This investigation was conducted at the United States Naval 
Postgraduate School, Monterey, California, during the period August, 
1964, through April, 1965. 

The authors wish to express their appreciation for the advice 
and assistance given them by Professor E. J. Andrews of the Aero- 
nautical Engineering Curriculum. Appreciation is also extended to 
the Ryan Aeronautical Company for making available the XV-5A data 


used in the analysis. 





2.0 DISCUSSION 


An investigation of the longitudinal and lateral-directional dynam- 
ic stability of the XV-5A VTOL aircraft in the transition flight regime 
was conducted by means of an analog simulation of the appropriate equa- 
tions of motion. The aircraft was taken to be in the transition configu- 
ration with 45° flap deflection and 15° aileron droop. Longitudinal po- 
sition of the center of gravity was taken to be at Station 246.0. The 
flight condition was assumed to be at sea level and out of ground effect. 

In the performance of a transition maneuver from hover to con- 
ventional flight, and vice versa, a fixed wing VTOL aircraft may en- 
counter the complete range of angle-of-attack -90°° <a < +90° due to 
the infinite number of possible combinations of vertical and forward 
speeds which might be achieved. In the present analog study of the dy- 
namic behavior of the XV-5A, the responses toa step or impulse per- 
turbation of pitch attitude and lateral velocity were measured for spe- 
cific steady state conditions of angle-of-attack and relative velocity. 
Time permitted a detailed investigation of only the positive angle-of- 
attack range 09 кос < 1909? corresponding to a descending transition 
from conventional flight to the hover. Some qualitative conclusions 
were drawn, however, concerning the dynamic behavior in the negative 
angle-of-attack range. 

Only the stick fixed dynamic stability of the XV-5A aircraft was 


investigated. Further, it was assumed that the symmetric and assym- 


metric motions of the aircraft were uncoupled, allowing separate 





analyses of the longitudinal and lateral-directional dynamics. 

Before the analog simulation could be conducted, it was neces- 
sary to obtain estimates of the aerodynamic characteristics of the air- 
craft over the complete angle-of-attack range -90° = x, = 1909. Reli- 
able aerodynamic data was available for the low angle of attack range, 
herein considered to extend from cx, =0° to cx 716^. (12, 13) For the 
range of large angles-of-attack, defined to extend from << =16% to ړم‎ 
=90°, the available low angle-of-attack data was extended by various 
methods described in detail in following sections. Estimates of some 
of the longitudinal stability derivatives were obtained from lift, drag, 
and pitching moment data resulting from static wind tunnel testing of a 
very simple model of the XV-5A aircraft. The wind tunnel tests of this 
model at large angles-of attack are discussed in Section 2.1. Estimates 
of the non-dimensional stability derivatives, presented in Section 2.2, 
were based ona system of stability axes, or wind axes, since the avail- 
able data for low angles-of-attack was expressed in this system. It was 
desirable, however, to use dimensional stability derivates referred to 
body axes in the analog simulation. The stick fixed, uncoupled equa- 
tions of motion, referred to body axes, are defined in Section 5.1. The 
relationships for converting the various non-dimensional derivatives to 
dimensional form and the equations for transforming these derivatives 
from stability axes to body axes are given in Section 5.3. Also includ- 
ed in this section is a digital computer FORTRAN program designed to 


perform these conversions. The analog circuits for the simulation of 
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the equations of motion are described in Section 5.4. The analysis of 
the results obtained is presented in Section 2. 3. 

Although the phugoid oscillation is usually assumed to occur at 
constant angle-of-attack in conventional stability analyses, it was not 
clear at the outset of this investigation that this would be the case for 
flight at angles-of-attack above the stall. Since most of the longitudi- 
nal stability derivatives vary with eX, in a highly nonlinear manner 
over some portions of the angle-of-attack range, it was considered 
necessary to simulate the variation of these coefficients by the use of 
variable base diode function generators in the analog computer circuit. 
The variations of the longitudinal dimensional derivatives for a devia- 
tion of angle-of-attack of + 1° froma given reference angle were pro- 


grammed into the function generators. 
2.1 Wind Tunnel Tests 


In order to obtain the basic lift, drag, and moment data of the 
XV-5A aircraft throughout the angle-of-attack range -90°= у= +90°, 
static wind tunnel “×2 were conducted on a scale model of the aircraft. 
The model, shown in Fig. 2.1, was constructed to a scale of one inch 
equal to two feet using data obtained from RAC. (12) A three view 
drawing of the XV-5A aircraft is shown in Fig. 2.2. No attempt was 
made to construct a sophisticated reproduction of the aircraft as far as 
airfoil sections, flaps, and controls were concerned. Rather the em- 


phasis was placed on obtaining the proper planform, scaled dimensions, 








and streamlined contours. It was anticipated that above the stall angle- 
of-attack the lift, drag, and pitching moment produced on the aircraft 
would be the result of a "flat plate' force effect since flow over all air- 
foil components would be completely separated. Planform would there- 
fore be most important. Below the stall angle-of-attack it was antici- 
pated that the data of RAC would be used in the analysis. (12) Never- 
theless, data was obtained in the region of low angles-of-attack in order 
to compare results with RAC as an indication of the adequacy of the 
model used. 

Model testing was conducted in the Aerolab Subsonic Wind Tun- 
nel located in Building 225 at the U. S. Naval Postgraduate School. A 
planform layout of the tunnel is shown in Fig. 2.3. A view of the tun- 
nel test set-upis given in Fig. 2.4 showing the complete model mounted 
by a three point suspension at a positive angle-of attack of approximate- 
ly 45°. Various sting mounts other than the one shown were used in the 
lower and negative angles-of-attack ranges. The mechanical balance 
used to measure the forces and moments is also shown. 

Model test data was obtained in both the WB VH and WBV con- 
figurations in the angle of attack ranges of o? <с<, < 4829 апа -829 < 
کے وپ ہش‎ -43°. Due to the physical constraints of the mounting system 
employed, angles of attack above 82° in magnitude were unattainable. 
Data in the range -43°=&,, = 0? was not obtained because of sched- 


uling difficulties of the tunnel facilities. 


Initial testing was conducted in the angle-of-attack range 








429 < کیو کچ‎ 82° at three flow velocities: 30, 60, and 120 feet per 
second. Preliminary analysis of the data obtained showed no notice- 
able variation of the results between the various flow velocities. The 
remaining tests were conducted, therefore, at 120 FPS. 
All tests were conducted in a clean aircraft configuration. 
Results of the tests are presented in non-dimensional form in 
Figs. 2.5 and 2.6. Fig. 2.5 shows Cy, and C, as functions of wing 


angle-of -attack, Xy, While Fig. 2.6 shows Cm.,, as a function of eX. 


Ë 
Both figures include RAC data for comparison. (12) 

In view of the relative crudeness of the model used in the tests, 
agreement with RAC data is considered good, particularly for Cmo 
The lift curve slope is in very good agreement but wing stall occurs at 


lower values of ex у and essc The Cp data 1s most in error since 


x 
no attempt was made to correct for flow into the engine intake duct. 
Consequently, the duct of the model presented a flat plate surface to 
the flow which produced a rather severe drag increase in the low angle- 
of-attack region. 

Above the stall angle-of-attack, both Cry, and Cp drop off 
abruptly indicating flow separation and a breakdown of the circulation 
lift mechanism. This decrease continues to an angle-of-attack of ap- 
proximately 20? when both Cr, and Cy begin to rise once again. How- 
ever, above the stall angle-of-attack lift and drag forces are predom- 


inately force components, along specified axes, of the total force on 


the aircraft resulting from a 'flat plate drag! effect. This is borne 





out by the fact that Cj, and Cp are equal at approximately < Q =+46° 
and Cj, reaches a maximum value, above the stall, at approximately 
ОХ „<+4710. A pure 'flat plate' effect would produce these effects at 
ex, 1459. To be noted also is the high value of Cj, reached in this region. 
EE value of Cr, equal to 0.9 at مه‎ ۳ is only approximately ten per- 
cent lower than C], moy at the stall where the wing is producing lift 'ef- 
ficiently'. As ex approaches ШОО; Cr, approaches zero indicating 

the flat plate nature of the results. 

In the angle-of-attack range +36°S eX $1409 Cpwypy is coin- 
cident with CLwpB vH indicating a complete blanking of the horizontal 
tail. On either side of this angle-of-attack range, the Cy yg, and 
CLWBV curves diverge. As ex approaches %909, however, CLwpvy 
approaches zero as did C{LWwpyy- The appearance of the deep stall 
region where the horizontal tail is completely ineffective is equally ob- 
vious in both the Cp and Cmcg plots. In this region of angle-of-attack 
1t may be expected that both the static and dynamic stability of the air- 
craft would deteriorate sharply. The static stability deterioration is 
most vividly shown in the Cmcg plot. 

Although data in the angle-of-attack region - AS? = ote. <- 09 
was not obtained, it can be expected that the deep stall region would 
not appear in this range since the horizontal tail would not be blanked 
at any negative angle-of-attack. This result is indicated by the avail- 
able data in that CLWBVH ‚ CDwpvp; and Cmcg WBVH are signifi- 


cantly different from Стуву, СОуву: and Стаср WBV in the region 





of = - 430 and show no tendency toward coincidence. Once again 
Cy, and Cp become equal at &,,= - 449 with CT, indicating a peaking 


- 439, 


tendency at approximately <<, 

Finally, above the stall angle-of-attack, it can be seen that the 
slope of the Спасе curve is such that a high degree of static stability is 
inherent in the aircraft throughout both the positive and negative angles- 
of-attack except in the deep stall region of positive angles-of-attack as 
previously discussed. 

With this data available it was possible to predict the static 
aerodynamic characteristics of the aircraft in the transition configura- 
tion throughout the angle-of-attack range -909 € ex,, S 4909. The final 
CLWBVH 2nd Cpwpvyg curves used in this analysis are presented in 
Fig. 2.7 while the final Спасе wBvp Curve used is shown in Fig. 2.8. 
These curves are made up, in general, of three distinct parts as indi- 
cated. For angles-of-attack 00£<<x, < +179, RAC data is used. (12) 
For angles-of-attack 4179€ e«.. € 182? and -839Sa«,, € -43? the wind 
tunnel data just discussed is used. For the remaining angle-of-attack 
regions, estimated variations have been used. These areas are pri- 
marily the negative angles-of-attack from 0% to -43% and that portion 
of the lift curve ТИЕТ following the positive stall angle-of-attack. 
Due to the 'flat plate' effect above the stall, the presence of flaps and 
drooped ailerons will have only a minor effect on the Cy, and Cp data 
obtained from the wind tunnel tests without these features incorporated 


on the model. For this reason, the data from these tests was used 





directly for the aircraft in the transition configuration. 

Fig. 2.9 presents the variation of Clt and Cpt, based on hori- 
zontal tail area, as a function of angle-of-attack. The variations for 
angles-of-attack 31702 се 2 +820 апа -839 < œ< Ж -430 were ob- 
tained from the differences (CL wpyp-CLwBv) and CDWBVH-CDwBv) 
of Fig. 2.5 and converted to a basis of tail area. These portions of the 
plots are presented as a function of wing angle-of-attack,™,, while 
that portion below the stall, obtained from RAC data, is presented as a 
function of the tail angle-of-attack, < . (12) Use of the data in the latter 
region for prediction of the effect on the aircraft as a whole requires the 
standard corrections for 9€/3x and 7 to establish the actual angle- 
of-attack and Cy, of the tail. Data above the stall angle-of-attack may be 
applied directly to the aircraft as a whole since tail blanking effects are 
obviously inherently included in the presentation. 

As in the case of the complete model, Cry and Cpt are equal at 
angles-of-attack of approximately +46% and - 43%, As << approaches 
90° Crt approaches zero while Cp; approaches a peak value. Above the 
stall angle-of-attack, Cr[, reaches a peak value at angles-of-attack of 
-430 and 4709. The tail blanking effect in the positive angle-of-attack 
range results in a lower peak value of Cy ¿ than in the negative angle-of - 
attack range and at an angle far removed from the expected 'flat plate ' 
Z5 


The composite curves of Fig. 2.9 are those used in the 


10 





following parts of the analysis. 
2.2 Non-Dimensional Stability Derivatives 


As a prerequisite to the dynamic stability analysis, it was 
necessary to estimate the magnitude of each of the relevant non- 
dimensional stability derivatives over the range of angle-of-attack 
from 09 to 909. Good estimates of all of the non-dimensional sta- 
bility derivatives were available for positive low angles-of-attack 
up to the stall at ex -16? (12, 13). Various methods, described in 
this section, were used to extend this data for large angle-of-attack. 
It should be noted that the lateral-directional rate derivatives are 
the most difficult to predict, either analytically or experimentally, 
even at low angles-of-attack. The complex empirical techniques 
used to estimate these derivatives below the stall are valid only in the 


range of angle-of-attack for which the slope of the lift curve, C 


Шад 
is constant (Z). In the absence of accurate means of predicting 
the values of the lateral directional rate derivatives at large angles- 


of attack, it was necessary to use elementary methods of analysis to: 


indicate their approximate variation over this range. 
2.2.1 Longitudinal Non-Dimensional Stability Derivatives 


II tati Longitudinal Derivatives, С ‚С ‚С 


L Dex 


Values of the derivatives с : e and for the WBHV and 
x 


m 


WBV configurations were obtained from the lift, drag, and moment 
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curves described in the previous section. These derivatives are 
plotted in Figs. 2.10 through 2. 12 inclusive for the positive angle- 


of-attack range, 0° < Сте £90%. Also included in Fig. 2.10 is 


E. - 


E» 1.2 C no 


La 


The lift and pitching moment coefficients due to rate of change 
of angle-of-attack arise from the lag in downwash at the horizontal 
tail. Equations (2-1) and (2-2), used in estimating یں‎ 4 T are 
derived in Section 5.2. 

en: 26 í. > 4 де (2-1) 
= e«t % 5, С dx 


ша 
Сава T (2-2) 


* 
۱ че 2€ . 

It has been shown that the slight variation of “3g with angle-of-attack 

is approximately linear below the stall for the XV-5A. (12) Since down- 

wash results from circulatory flow about the wing and since the wind 

tunnel tests, discussed in the previous section, showed that circulation 


decreased rapidly above the stall, it was concluded that C and 


L < 
С. e approached zero as angle-of-attack increased above the stall. The 
ہی‎ 
data of Ryan for C, ¿and C „were therefore extended to reach zero 
X ос 
at EC 289 and were replotted in Figs. 2.13 and 2.14. (13) Both of 


these derivatives were assumed to be qualitatively similar at positive 


and negative angles-of-attack. 
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2 ПОЕ ӨРЕ 
a q 

The wing and horizontal tail make significant contributions 

to the derivatives С, аала С: . Curvature of the flow relative to a 
+ 

translating wing occurs when a pitch rate, q , 18 imposed on the 
wing. This effect can be thought of as an equivalent change in cam- 
ber which alters the lift and pitching moment characteristics of the 
wing. The empirical relation of Equation (2-3) was used in estim- 


ating this effect (6). 


2а x E 
a,” (3 - 37 


The quantity xc varies over the angle-of-attack range according 


Ое geometric construction of Fig. 2.15. 


lhe contribution of the horizontal tail to the derivative results 
from the change in angle-of-attack due to the vertical velocity com- 


ponent of the tail and is given by Equation (2-4) which is derived in 


Section 5.2. 


C = РС 1 St 4 (2-4) 
The variation of Же with =< is shown in Fig. 2.16. The sum of the 
components C, and C were plotted in Fig. 2.17 over the positive 
Чу 9+ 
and negative range of angle-of-attack. For negative angles-of-attack 
the horizontal tail does not pass through a region of deep stall or 


separated flow from the wing. The dotted portion of the curve of 


Fig. 2.17 represents an estimate of C. in the range of negative 


angle-of-attack corresponding to the positive angle deep stall range. 
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The contribution of the wing to the pitching moment coefficient 
due to pitch rate, С» , was estimated by means of the empirical Equa- 
3 


tion (2-5). (6) 


А+е Сов Л. 2 A+PCo 
nq „В-< E x a بی‎ ta “Ле i 


The inboard and outboard sections of the wing of the XV-5A have 
different sweep angles. A mean sweep angle, A= 209, was assumed 
for use in the above equation. 

Ihe effect of the horizontal tail on C was estimated from 
Equation (2-6). 

C ma, = С, A (2-6) 

The derivative plotted in Fig. 2.18 against angle-of-attack. 
As in the case of "id m not altered by the deep stall condition 


in the negative angle-of-attack range, as indicated by the dotted 


portion of the curve. 
аа Се С И 
ли Б 


The non-dimensional derivatives of forces normal and paral- 
lel to the wind axis and of the pitching moment due to airspeed are 
given by Equations (2-7), (2-8) and (2-9) which are derived in Section 


5. 2. 


ЖЕ узу а (2-7) 
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AEB (2-8) 


Camu = ЕС. — 


Since the aircraft is assumed in trimrmed flight with thrust pertur- 
bations ignored, C was neglected. Also neglected in the deriva- 
u 


tion of Equations (2-7, (2-8) and (2-9) were Mach number effects 


which are of no interest in this study. 


2.2.2 Lateral-Directional Non-dimensional Stability Derivatives 


ФЕР‏ ادف د 


The sideforce due to sideslip results from imposing a lat- 
eral angle-of-attack, В , on the fuselage and on the vertical tail. 
The wing was assumed to provide a negligible contribution to this 
derivative. At low angles-of-attack a large resultant negative side- 
force, acting close to the nose of the fuselage, can be expected when 
a sideslip angle is established. As €x is increased, the effective 
sideslip decreases, the center of pressure of the sideforce will 
move aft on the body, and the magnitude of the sideforce will decrease. 
In the limit, as cx approaches 909, only a small sideforce, acting at 
approximately the geometric centre of the body, will be exerted. For 
angles-of-attack below the stall, Cy was obtained from experimental- 
ly determined curves. (12) This data was extropolated to an estim- 
ated small value of E at Es 90° as shown in Fig 2.19. 


Pp 
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A negative sideforce is exerted on the vertical tail by a side 


with €X, indicated in the experimen- 
۷ 


slip. A slight variation of Cu, 
tal results for low, is probably caused by fuselage interference. 
(12) As angle-of-attack increases to the stall, part of the vertical 
tail becomes immersed in the separated flow behind the wing. The 
effective tail area and C are decreased. This deterioration be- 

6v 
comes more severe as E further increased and more of the tail 
area is covered by the separated flow region. Since it appears that 
at least part of the vertical] tail will be in unseparated flow at all 
times, it was estimated that C would de creage to a minimum of 

6۷ 

15 percent of its average value for low angles-of-attack. It was 
noted that the vertical tail would have minimum effectiveness at a 
smaller angle-of-attack than would the horizontal tail. For this for- 
mer, this angle was estimated to be 329. As the vertical tail em- 
erges from beneath the wing, at very high angles-of-attack, it tends 
to regain effectiveness. However, for very large wing angle-of- 
attack, the effective lateral angle-of-attack produced by a given side- 
slip velocity will be very much greater than the stall angle of the 


vertical tail. Therefore C will be very small for large e 


da y 


Cu, at E c 459 was therefore estimated to be about 35 percent of 
۷ 


its magnitude below the stall. Above CX - 459, C decreased 


w de, 
toward zero at €« - 909, 
W 


The estimated variation of С, with CX is plotted in Fig. 2.109, 
E w 


based on available data for C at low anglé-of-attack. (12) 


4e 
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280.202 C ng 


The yawing moment due to sideslip comprises contributions 
from the fuselage and vertical tail. These moments are the products 
of the sideforces described above acting on lever arms about the air- 
craft centre of gravity. The moment arm of the fuselage sideforce, 
measured trom the fuselage center Of pressure aft to the centre of 
gravity, decreases to zero as یرت‎ 600 hes 909. The moment arm 
of the vertical tail sideforce is measured from the aerodynamic cen- 
ter of the mean aerodynamic chord of the vertical tail. Due to the 
location of this point relative to the aircraft center of gravity, the 
vertical tail moment arm will not vanish at p 909, but will be re- 
duced to a minimum of about 20 percent of the arm at کے‎ 0°. 

The wing contribution to с. was considered to be negligible. 

Fig. 2.20 shows the extrapolation of the variation of C. at low 


ё 


angle-of-attack to the large angle range (12). 


282110053 Cig 


The dihedral effect of the wing and the sideforce acting on the 
vertical tail comprise the rolling moment due to sideslip. Experi- 
mental results show that C4 increased negatively with angle-of- 

w 
attack up to the stall, apparently due to fuselage interference. (12) 
Above Te stall, an abrupt reduction, in differential lift on the wings 


would be expected for a given sideslip. The reduction in lift results 


from the decrease in effective angle-of-attack change due to dihedral. 
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The ampoent шу, > is the productief C and the moment 
V ۷ 
arm measured perpendicular to the relative wind axis, through the 
center of gravity, from this axis to the aerodynamic center of the 
vertical tail. For small angles-of-attack the aerodynamic center of 
the tail is above the relative wind axis and a negative rolling mo- 
ment results from a positive sideslip. For large positive, 
however, the position of the vertical tail is such that a positive 
rolling moment is caused by a positive sideslip. The length of the 
moment arm and therefore the value of Cle are zero at x = SEN 
V 

In the гапре 169< X ,= 36°, Clo, remains small due to blanking of 
the vertical tail by the wing. For very large = although the side- 
force on the vertical tail due to sideslip would be small, a signifi- 
cant rolling moment would be expected due to the very large mo- 
ment arm on which the sideforce acts. Note that the discussion 
refers to stability axes and that a rolling motion, about the X sta- 
bility axis at < 90° has become a yawing motion about the 2 
body axis. 

For the negative range o the moment arm will always be 
such that mo be negative. No blanking of the vertical tail would 
be experienced in the range -169< ex = - 36° as in the positive range. 
Rather, it would be expected that at very large negative angles-of- 
attack, the vertical tail weguld become progressively blanked by the 
horizontal tail. 


The available data for Ey at low angle-of-attack was extended 
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to large angles by means of the above methods and the resulting 


curve is presented in Fig. 2.21 (12). 
pus. 2.4 
Ca, 


The major contribution to the sideforce due to yaw rate, Cy 
comes from the vertical tail. The empirical methods mentioned 
previously were used to estimate the very small contribution of the 
wing to the derivative Cy at low angles-of-attack. (13) The effects 
of wing and fuselage on Cy و‎ assumed to be negligible at large 


angles-of-attack. For the range 90 < 90°, Equation (2-10) 


was used to estimate Cy . (7) 
Ary 


C = P C (2-10) 


dry "By 


Fig. 2.22 shows the variation of Cy with Oy. 
۲ 


SO, 


The yaw damping derivative, Cy, results from differential 
drag due to differential velocity over the wings and from the yawing 
moment due to the lateral angle-of-attack change associated with the 
yaw rate on the vertical tail. Equations (2-11) and (2-12), derived 
in Section 5.2 for assumed constant span wise lift and drag distri- 
bution, were used to estimate the wing and tail contributions to 


at large angles-of-attack 


С. = <> (0.8) (2-11) 
3 


LN 
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el © (Dale 
Спа E "p, 


The variation of Сы, with ex, is shown in Fig. 2.23. 


2.۰2 6 C e, 


The rolling moment due to yaw rate generated by the wing 
results from differential lift due to velocity changes over the wings. 
Handbook methods were used by the manufacturer to account for 
geometric effects and non-constant span wise lift distribution in estim- 


ating Со, for small angles-of-attack. (2,13) For large с 
W 





W 2 
Equation (2-13), an approximate relationship, was used. 
S 
E = —= (0.8 (eau) 
به‎ 2 ( ) 


The derivation of Equation (2-13) is given in Section 5.2. 

Because of lack of circulation at large angles-of-attack, the 
assumption of constant lift distribution used in deriving Equation (2-13) 
was considered to be more accurate than below the stall where circula- 
5977٠6 ٦16 15 21۰ء6‎ 

The vertical tail contribution to Cy was computed from 


Equation (2-14). (7) 


@ = E Саа, Copy (2-1) 
PON C Jg y 


(o is plotted against Ow for positive and negative angles in 


Fig. 2.24. 
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DNDC2.7 C yp 


The only significant component of the sideforce derivative 
due to rolling moment comes from the vertical tail. Equation (2-15) 


ine Cun. 
was used to determine 1e (7) 


OG (2-15) 
Coe = ma + C Ее 
Яр tey © “Яв, ӘР 
29; 
The second term of Equation (2-15) describes the effect on the verti- 
cal tail of sidewash due to roll rate. As angle-of-attack was in- 
ه‎ 
creased above the stall, age was assumed to approach zero rapid- 
U 


o 
ly. At low& „there exists a small sideforce due to roll rate which 


arises from wing dihedral. This effect is very small compared to 
the vertical tail contribution and since it decreased with xit was 
neglected. с is plotted in Fig. 2.25 for the range - 909 € =< N = 


+ 
900. 


2 RE Cas 


Both the wing and the vertical tail have important effects on 
the yawing moment due to roll rate. The contribution of the wing 
comes from the change in magnitude and inclination of the resultant 
force vector on each wing. From elementary analysis considering 
constant span wise life and drag distribution, Equation (2-16) was 


obtained. 


Chp, = (Ca. = м ES m 


Derivation of this expression is shown in Section 5.4. For very low 
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angles-of-attack, this elementary analysis is in reasonable agree- 
ment with experimental data. (13) At very high angles-of-attack, 
where circulation no longer exists, the above expression is also 


considered valid. It is important to note that C must be obtained 
A 


experimentally for high angles-of-attack since induced drag results 


from circulatory lift only. Therefore results of wind tunnel tests 
previously discussed were used in evaluation of ы high angles- 
of-attack. The method just described predicts negative values of 
for angles-of-attack just above the stall where some circulation still 


exists. Experimental results for wings of similar planform indicate 


that Ch, will be positive in this region in contradiction of the theory. 
W 


(8) Therefore from the data of this reference, Chp was estimated 
іот 1605 с<,, < 329, This estimate was faired into the very high 
angle-of-attack estimate and into the low angle data previously 
mentioned. 

The vertical tail contribution was determined from Equation 
Z217). -(7) 


C C Әс 
2D DEN Ay Ip NE. (2-17) 
Ce, C 48 y ۷ oF 


The variation of Chpwith = „is shown in Fig. 2.26 for - 90%<—¿ < 909, 
px cC 
Ap 


Differential lift on the wings due to differential change of 
angle-of-attack provides the largest component of the roll damping 
derivative. ha estimated from Equation (2-18) for large 


angles-of-attack. 
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ch С). (23 te) 


The derivation of Equation (2-18) is given in Section 5.2. The contri- 


bution of the vertical tail, given by Equation (2-19), was found to be 


negligible. (7) 


C (ley oc 

D + Cy —‏ ے 

C O>. Счву Ав, ot 

was plotted in Fig. 2.27 for the range -90° < O SR 





(2-19) 


2.3 Analysis of Analog Computer Results 


The longitudinal response of the aircraft to perturbations in 
lateral velocity were measured for the discrete reference flight con- 
ditions of steady state angles-of-attack and relative airspeeds listed 
in Tables IX and XIII in Section 5.4. All computer runs were made 
for a steady state pitch attitude of O, = 10° except for one run at 
Е 80, V5 2100 knots, and 5) ج‎ 09%. In the longitudinal analysis the 
perturbations U, n сх, Ө, апа q were recorded on a multi-channel 
strip recorder. Perturbations of Y”, dp, Y, and ¥ were recorded 
in the lateral-directional tests. 

The first computer run of each of the longitudinal and lateral- 
directional series of tests was conducted for eT 8%) Vo= 100 knots, 
and e 109. These reference conditions correspond approximately 
to the minimum speed for which the period and damping parameters 
obtained from these initial computer runs were compared with the data 


of RAC for trimmed conventional flight at V5- 94. 6 knots and e - 8. e, 
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For the longitudinal case, from Fig. 5.12, the period of the 
phugoid oscillation was measured to be 39.5 seconds, exactly the 
period given by the manufacturers report. (13) The measured time 
to damp to half amplitude was 20.0 seconds for the initial run at 
o. = 10%, A second run with Ə = 09, but otherwise similar to the 
first, resulted in a value of T 1/2 = 12.3 seconds measured from 
pig. 5.15. These values of T j/2 are in reasonable agreement with 


the contractor value of T ]/2 = 14.6 seconds for ӨЯ - ٦ 

Similarly, the initial lateral-directional test produced a 
period and time to damp to half amplitude for the dutch roll oscil- 
lation of 3. 75 seconds and 2.5 seconds respectively, as measured 
{rom Fig. 5.35, which compare very favorably with the КАС дата 
of 3. 32 and 2. 29 seconds respectively. (13) In addition, the dutch 
roll parameter,|$% |, was measured to be . 612 compared with the 
given value of . 610. 

Once the analog circuits had been proven in the manner 
just described, the computer runs were carried out for the refer- 
ence flight conditions of interest. The results of each run are 
shown in Figs. 5.12 through 5. 34 for the longitudinal mode and 


EES. 5. oo through 5.5/ tor the lateral directional mode. 
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2.3.1 Longitudinal Mode 


One of the most significant conclusions which can be 
drawn immediately from an examination of the longitudinal results 
is that the deviation of the angle-of-attack from its reference 
flight condition is very small for all runs. The variation is so 
slight that it is concluded that the longitudinal stability derivatives 
could be considered constant at any given reference angle-of-attack. 
Hence, in a future simulation of this nature, the function generators 
could be omitted from the analog circuit. The stability derivatives 
would be inserted as simple potentiometer settings thus simplifying 
the circuit. 

It was found that for angles-of-attack at which the longi- 
tudinal response was oscillatory, the period varies considerably as 
a function of both angle-of-attack and airspeed. Table I lists the 
period and frequency of the longitudinal response at the various 
reference flight conditions. This data is not presented graphically. 
Such a presentation would be misleading since ranges of angle-of- 
attack exist where the motion is aperiodic. From Table I, however, 
it is clear that the period becomes quite short in the vicinity of the 
stall, increasing at very high and low angles-of-attack. It may also 
be seen in Table I that at low angles-of-attack the period decreases 
slightly with reduced airspeed while at large angles the opposite is 


true. 
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Wide variations in the damping characteristics were ob- 
served over the ranges of angle-of-attack and airspeed. For the 
purpose of graphical presentation of these variations, a damping 
parameter l/T was defined. The quantity, T, refers to the time to 
half amplitude, in the case of a stable oscillation, or time to double 
amplitude, in the case of an unstable oscillation or pure divergence. 
B T2 was measured for each flight condition from Figs. 5.12 
through 5.34. The parameter 1/T plotted as a function of angle-of- 
attack is presented in Fig. 2.28. The value of 1/T = 0 in this figure 
represents the condition of neutral stability or undamped oscillation. 

At low angles-of-attack it is clear from Fig. 2.28 that, 
as airspeed decreases, the damping, and hence stability, deterio- 
rates. At P s 89 the aircraft is stable at 100 knots, neutrally sta- 
ble at 60 knots, and unstable at 20 knots. At the stall there is a 
sharp increase in the rate of divergence for all airspeeds. It is 
important to note that, in contrast to the low angle-of-attack region, 
the rate of divergence is greater at 100 knots than at 20 knots for the 
range of ST above the stall. This phenomenon occurs because of 
the change in polarity of major longitudinal derivatives such as MET 
and Mg near the stall. It is not possible to be more specific about 
this point because of the different manner in which these derivatives 
vary with angle-of-attack. 

The second region of rapid deterioration of stability, shown 


by Fie 2.268, iS ih the vicinity of E 289. This deterioration results 
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from a combination of several factors. At the stall the derivative 
Mg becomes positive and remains so up to an angle-of-attack of 
approximately 349. The derivative M remains negative until an 
angle-of-attack of approximately 269. It then becomes positive and 
combines with Ma causing the rapid divergence indicated at = 269. 
At an angle-of-attack of approximately 40° М again becomes nega- 
tive thereby improving the longitudinal response considerably as 
seen in Fig. 2.28. The fluctuations of Mq and M ‚just described are 
the direct result of horizontal tail stall and subsequent blanking in 
this region of angle-of-attack as discussed in Section 2.1. 

The longitudinal motion of the aircraft is slowly divergent 
at very large angles-of-attack for the practical range of airspeeds. 
The trend of stability E airspeed variation is reversed once again 
in that the instability becomes more severe as airspeed is decreased. 

In general, it is concluded that at low and very large angles- 
of-attack for all practical airspeeds of the transition flight regime the 
rate of divergence is slow enough that a pilot could maintain longi- 
tudinal control of the aircraft. However, in the range 16% SX 838" 
the rate of divergence is such that the pilot could maintain control 
only at a reduced airspeed. If a completely flexible ind angle- 
of-attack capability is required or desired an automatic stabilization 
system would be necessary in the longitudinal control loop. 

Because of the limitation of time, an analysis of the sta- 


bility in the negative angle-of-attack range was not carried out. 
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However, some qualitative conclusions can be drawn about stability 
in this range from the observations of stability in the positive angle- 
of-attack range. At low and very large negative angles-of-attack 
the longitudinal stability characteristics of the aircraft could be 
expected to be similar to those already discussed for positive angles- 
of-attack. The stall would occur at a smaller negative angle than 
that positive angle at which it occurs but this would not alter the 
nature of the dynamic response of the aircraft appreciably. The 
rapid deterioration of stability in the deep stall region of positive 
angle-of-attack would not occur in the negative range since the hori- 
zontal tail would never be blanked by separated flow from the wing. 

For a number of reference flight conditions for which it 
was felt that a significant effect might be noticed, the nose fan pitch 
damping was removed from the analog circuit. Fig. 2.28 shows 
that a very slight reduction in damping resulted. There was notice- 
able change in period of oscillation with the removal of nose fan 
damping. 

Main fan lift damping, although included in the analog 
circuit, was known to be small in comparison to the wing damping 
and would not change the basic characteristics of the aircraft. 
Therefore computer runs were not performed with fan lift damping 


removed. 


2.3.2  Lateral-Directional Mode 
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The periods and frequencies of oscillations of the lateral- 
directional responses are presented in Table II as determined from 
the analog recordings of Figs. 5.35 through 5.57. The period of the 
lateral-directional oscillatory response varies directly with angle- 
of-attack and inversely with airspeed. This dutch roll motion ex- 
hibited periods ranging from 3. 75 seconds at V¿ = 100 knots and 
هت‎ 9» to 20 seconds at Vo - 20 knots at ex, - 729. Considering 
the periods of the longitudinal responses and those of the lateral- 
directional mode above, it can be seen that all periods are consider- 
ably larger than one second. It is concluded, then, that pilot induced 
oscillations caused by coupling of control deflections with the natural 
frequency of the oscillations would not constitute a problem in the 
transition flight regime for this aircraft. 

The lateral-directional damping characteristics of the air- 
craft over the ranges of speed and angle-of-attack are shown in Fig. 
2.29 by means of the 1/T parameter as described in the previous 
section. In this case it is clear that there is a deterioration of sta- 
bility with reduced speed in most of the angle-of-attack range. Ex- 
cept for a small range of angles-of-attack just above the stall, the 
degree of this deterioration is not as severe as it is in the longitud- 
inal mode. It is noted that, generally, the aircraft will not become 
unstable in the lateral-directional mode until speed is reduced to the 
vicinity of 40 knots. 


Within the small range of angle-of-attack, 16° to 19°, the 
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lateral-directional motion becomes rapidly divergent because of the 
change in polarity of the important roll damping derivative Cap: 
Comparison of Figs. 2.28 and 2.29 indicates that the peak rate of 
lateral-directional divergence for V = 60 knots is approximately 
double that of the longitudinal divergence at the same speed. 

An interesting conclusion drawn from Fig. 2.29 is that 
blanking of the vertical tail in the separated flow behind the stalled 
wing, that is in the angle-of-attack range 24° © s 409, seems to 
have a serious destabilizing effect only at higher speeds. At 20 knots 
blanking of the vertical tail appears to have little or no effect on the 
lateral-directional response. 

Conclusions similar to the longitudinal mode are drawn con- 
cerning the negative angle-of-attack region as far as the lateral-dir- 
ectional motion is concerned. Blanking of the vertical tail would not 
occur and hence no deterioration of stability would be experienced due 
to this effect. However, the rapid divergence in a narrow angle-of- 
attack band immediately following the stall in the negative angle region 
could still be expected to occur. 

Main fan roll and yaw damping were removed from the ana- 
log circuit for some runs for which a significant effect might be ex- 
pected. In only one case, at ox = 28? and V, - 60 knots, did the fan 
damping prove to be important. Under these conditions the time to 
damp to half amplitude of the dutch roll oscillation was increased 


by a factor of three. 
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1 158 evident irom the analog traces of Figs. 5.35 through 
5.57 that the aircraft exhibits a very slow spiral divergence. This 
minor instability was not rapid enough to be of interest. 

From the observed lateral-directional stability charac- 
teristics it is concluded that for flight.at any transition speed over 
the positive and negative angles-of-attack ranges, except +160 < 
A, se 209, (һе rates of divergence are such that a pilot could 
maintain lateral-directional control with a manual control sytem. 
Within the range +16%< eo ae the very rapid rate of divergence 
of the dutch roll oscillation would prevent a pilot from maintaining 
control with A manual control system. Therefore, sustained flight 
in this range of angle-of-attack would make necessary automatic y 
stabilization of the lateral-directional control loops. It.is suggested, 
however, that a pilot could transit through this critical range with | 
manual controls provided he remained within this angle-of-attack 
range for only a short period oí time. 

The dutch roll parameter, |, |, which indicates the rela- 
tive magnitude of the roll motion to the yaw motion, | was measured 
directly from the analog traces. The results are plotted as a func- 
tion of angle-of-attack in Fig. 2. 30. At low angles -of-attack | V | 
increases as speed decreases due to the deterioration of roll damp- 
ing with reduced airspeed. For angles just above the stall Б 
becomes much greater than unity. Noting, from Fig. 2.29, that 


the aircraft is very unstable in this range of angle-of-attack, it is 
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concluded that the instability is primarily a roll divergence. For 
large angles-of-attack, Fig. 2. 30 shows nat] % [increases with air- 
e 

speed. At 60 knots in this angle range, where the aircraft is stable, 
the rolling motion will be equal to or greater than the yawing motion 
whereas at 20 knots, where the aircraft is unstable, the yawing mo- 
tion will be more pronounced than the rolling motion. 

As a final indication of the regions of 0ء9‎ flight 
which would be hazardous with manual controls only, Fig. 2.31 is 
presented. This figure, in which relative airspeed is plotted against 
angle-of-attack with time to double amplitude, T>; as a parameter, 
is a cross plot of the curves of Figs. 2.28 and 2.29. Sustained 
flight within the speed-angle-of-attack regions outlined by the dashed 
boundaries is considered undesirable if manual controls alone,are used. 
The particular boundary used to define this prohibitory flight regime 
is dependent upon the minimum value of Time to Double Amplitude, 
T5, with which a pilot can operate. It can be seen that the oe 
directional boundary is very narrow in comparison to the longitudinal 
boundary although the rate of divergence is faster within the former, 
as previously discussed. The boundaries are approximate in nature 


and are not meant to indicate exact areas. 
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3.0 CONCLUSIONS 


1. Lift and drag forces at large angles of attack result from 'flat 
plate' type forces with no contribution from circulation effects, 

2. Blanking of the horizontal tail causes rapid deterioration of static 
and dynamic longitudinal stability in the range + 20%=cx, =+ 400, 
Blanking of the vertical tail causes a moderate deterioration of lateral- 
directional dynamic stability in the range + 24^ 5 eX S * 40? at the 
higher transition speeds. 

3. Longitudinal stability derivatives can be considered constant for 
any given reference angle-of-attack. 

4. Periods of longitudinal and lateral-directional oscillations over 

the ranges of speed and angle-of-attack considered are, in all cases, 
greater than 2. 5 seconds. Therefore, the danger of pilot induced 
oscillations in the transition flight regime is not considered to be a 
problem for the XV -5A. 

5. At low and very large angles-of-attack there is a moderate decrease 
in longitudinal stability with i$ em reduction. In the range + 160 > 

yw =+ 36° the aircraft exhibits a severe longitudinal instability, the 
rate of divergence varying directly with airspeed. The shortest ob- 
served time to double amplitude was 0.5 seconds at - +28% and Vo= 
100 knots. 

6. At low and very large angles-of-attack there is a moderate decrease 


in lateral-directional stability with reduction of airspeed. In the range 
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шо == = 124? the aircraft exhibits a very rapid roll divergence. 
The shortest observed time to double amplitude was 0.45 seconds at 
er * 189 and V, = 60 knots. 

7. Except for the range of angle-of-attack + 16°S x, = + 409 rates 
of divergence in regions of instability are such that a pilot could 
maintain manual control about all three axes. For sustained flight 
in the excepted range, automatic stabilization of the aircraft in all 
control loops would be required. 

8. In general, the contributions of nose fan pitch damping and main 
fan roll and yaw damping to the overall damping characteristics of 
the aircraft are small. 

9. The variation of the Dutch Roll Parameter with angle-of-attack 
indicates that the lateral-directional instability in the range + 169 € 
C. S4 249 is primarily one of roll divergence. The maximum 


observed value of | Pre | was 2.32 at Е - + 180 and Vo = 60 knots. 
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TABLE I 
Periods and Frequencies 


Longitudinal Mode 


Period C2, 
Sec. Rad/Sec 
С 0161 


39.0 










Dashes signity no oscillaticn recorded 
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TABLE Tr 
Periods and Frequencies 


Lateral = Directional -Mode 


Period 
ЩЕ сес 


LD4 









.394 
999 


Dashes signify no oscillation recorded 
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- Planform Layout - Subsonic Wind Tunnel 


Figure 2.3 
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SO APPENDIX 


5.1 Aircraft Equations of Motion 


The equations of motion used to define the dynamic behavior 
of the aircraft after being preturbed from a reference flight condition 
are those commonly used in a dynamic stability analysis. The devel- 
opment of these equations is presented in detail in many references 
and will not be shown here. (3,9,10) Only the final equations and 
the assumptions upon which they are based are presented. The sym- 
bology used is that of Etkin. (3) The analysis considers only stick 
fixed conditions; therefore no contributions arising from control 
deflections are included. 

The equations of motion were developed utilizing an Eulerian 
set of axes fixed to, and moving with, the aircraft mass center. For 
ease of definition of motions along and about these axes, a body axis 
system, as defined by Etkin, was chosen. (3) Fig. 5.1 presents the 
axis system and defines positive directions of forces, moments, and 
motions. The assumptions made in defining the equations are: 

l. The aircraft mass is considered constant. 

2. The aircraft is considered a rigid body. No aeroelastic 
effects are included. 

3. QGyroscopic effects of spinning masses are neglected. 

4. Higher order acceleration terms are neglected. 


5. The motion of the aircraft consists of small deviations from 
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some reference flight condition. 

6. The disturbance quantities and their derivatives are small 
so that squares and products are negligible in comparison to first 
order quantities. 

1. The ^ 7, plane is a plane of symmetry. Products of inertia 
involving the% y and 44 directions are therefore identically equal to 
Zero, 

The reference flight condition chosen was that of symmetric 
steady state motion. This signifies the absence of any steady state 
linear acceleration and rotational velocity or acceleration. Also 
inherent in this reference flight condition is the fact that Vo, ф. ; 
and Y are identically equal to zero. 

Based on the aforementioned conditions and assumptions, the 
six equations of motion may be divided into two sets of equations 
defining motions which are independent of each other. These equa- 


tions may be expressed as follows. 


Longitudinal Mode 


И 


A X - (md cos ره‎ т (а + weq) (S21) 


m (er - 4.4)  (5.2] 


il 


AZ - (m34 Sin ©.) Ө 


AM = B4 (5-3) 


Lateral-Directional Mode 
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AY + m q ( $ Cos O, + Y sin &)- m(v- + UN — Uk, 2 (9093, 
AL =A p -Er (5-5) 
4٣۸ Ep ) 


The force and moment terms AX, AY ; AZ; AL,AM, and 
AN consist of a sum of contributions from various aerodynamic 
effects only. 

Generally the above equations are completely non-dimension- 
alized permitting the previously mentioned forces and moments to be 
expressed in terms of non-dimensional stability derivatives. How- 
ever, this analysis was to be conducted with the aid of an analog 
computer using dimensional derivatives throughout the simulation. 
Hence, dimensional derivatives were defined as shown in Section 5. 3 


and the final equations of motion used are as follows. 
Longitudinal Mode 
X, v +X. u- -(mq соѕӨ,)Ө= т( + их, q) 20 
ZU + Zur 2 4244 + 4 ~0 = (i - 1.4) (5-8) 
Ми حثا ث٢ + جا‎ + G + АМ = B q (SS) 
Lateral-Directional Mode 


ул +Ьр +۷۱۴۱ mq( $cos &, « y sine ye m(v «us r- ur p) (5-10) 
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u» 





L. Ep «LU р - Ар - Ег (5-11) 


Ny + Мр + Маг + Nr r GED (5-12) 


As a final simplification, engine and fan RPM was taken as 
being invariant during a perturbation in forward or vertical speed. 
The fan damping terms arise primarily from a change in air mass 
flow through the fans as the aircraít is perturbed from its reference 


flight condition. 
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5.2 Simplified Relations for some Non-dimensional Stability 


Derivatives 


In Section 2.2 it is stated that some of the estimates of stability 


derivatives were based on simplified analyses. 


The derivation of 


expressions for these derivatives, obtained from an unpublished ref- 


erence, are shown in this section. (1) 


5.2.1 ам 
+ 


Standard expressions for ته‎ and € may be written: 


Soa - X 


€ = (xw = 42 at) 


The quantity At is the time required for the flow to travel from the 


wing to the tail and may be expressed: 


m 

At U 
>e | 
Тһеп: x= < li- 26 \ + dE de lt 
K t F е) 3x dt لا‎ 


From the above expression it can be seen that AX due to downwash 


lag is: 


A په‎ 


i 
ТТ» 

fh 

o 


S de kh ۔‎ ۷٦ 
Ü Oe C 


و 


A L 


5/2 075, (+ Ci 


AC. 7 £C, n, 2t 2€ (+ oz 
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EXC _ АС + 7 C S Р дє 

Finally: و‎ 3 = ES 5-13 

шашу: 64, 573 E (5-13) 
20 

Ze Ge 


EE 
This derivative results directly from the product of the left 


on the horizontal tail due to downwash and the tail arm. 


а a ^ 
ROMS E (5-14) 


pU с, 
dt 


The derivativeC,_ represents the lift on the horizontal tail 
caused by the change in angle-of-attack at the tail due to a pitch 


rate. The development is as follows. 


PER 7 


+ 
AC, = pU St C 21 C 
سک ا مع‎ + tor “E FU 


Finally: C, = (5-15) 


5 14 Cx, 


Referring to Fig. 5.2, the following expression may be 


written. 


76 





Assuming | invariant with u and differentiating with respect to خا‎ 





gives: 
B. 2 a) 
Su 7 ا 2 اد‎ 
E ar U 
C y T 
Cyu » š d ч u. [x D Su 
(%) A) 2/2 
шегеге: C. - —2 Cg (5-16) 
5 С 
ди 
Cain referring to Fig, 522; the КОО ae development can 
be made. 
Z =-Tsına —L 
22 SES “03 +8٣٦ 
0ئ‎ 6 ст z Е 


| 
| 
x 
> 
С 
(> 


С. 2-3 Ca n E ва a UA 
^" — M ч) nf > | 
Then: 8 = - 2 C, (5-17) 





5: 2.6 
С, 
From Fig. 5.3, the following relation may be written: 
AD = C, L (U-ry) cd 
Dw 2» d 


Then: 65 = 5‫ۀ‫ گ‎ Co, pu "S Ay 


Assuming a constant spanwise value of co gives: 


3 (AN) 2 = Co, / Оц с Ay 


r 
3 +b, 


T vf C4'd 
m bw /^ by d 4 


. 5 E 1 X (^ E 


Mo AU] 
3r CL 42V Sb b 
b. 
С U cu a Qu 
وت‎ E 4 24 A PUT SL 


cb‏ + مه ہ+۔ ہہ ле‏ پ9 
س — 


E 


= гама и 
S b> Sb 
Finally: kx. - - Ср, 0.8 (5-18) 
w 3 


The factor 0.8 improves the approximation for the area moment of 


inertia, Ixy. 
ва Е ~C 


The vertical tail contribution can be expressed as: 
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C C и 


е ٤ وپ‎ 20 


کم 0,۷ م » 


Fig. 5.4 defines the vertical tail moment arm. The angle of sideslip 


at the vertical tail may be expressed: 


A = ry 2 2, 





2% . rb 
U b 2.0 
Тһеп: Cno, E Eng 21, (5-19) 
V 
5.2.8 C 
Í 


Referring once again to Fig. 5.3, the following relations may 


be written: 


аі = с, ^ /(U- rqY eds. 


жр - Сър (- 204) сац 
7 


XL) = ا ےپ ا‎ 
> Co C, e U ç c dg 


Assuming constant spanwise lift distribution gives: 





i by " 
ыы се]. PM 
-b, 
er = C, p U N 
m с дере T NONU 
H 2 ya ms طک لا ما‎ b 


23 





& = +С, Ду 
hie, Sb 


Finally: Cy ME 6 (5-20) 
rw 3 
The factor 0.8 is again used to improve the area moment of inertia 


estimation. 


EE C 
N Pw 
Fig. 5.5 presents a spanwise element of the right wing under- 


going a positive roll rate. From this figure, the following develop- 


ment is possible. 


dN = yÍ (D+aD) cos ac - D- (L+Al)sına «| 


$ 2D=Lax) 


4 CCS, - с \ зоо е dy 40 


+b,, 
n 1 b 
N = ۲ XN |, 43 +6 2 
Cy = UC a), Da $e £5 
C = X Ca = 2 GG ~ С ТЕЕ 
ШЕКТЕ? єч 
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Finally: Cno = (Co, сө (5-21) 


6.0 


5.2.10 Су, 
یہ‎ 


From Fig. 5.5, the wing contribution to roll damping may also 


be analyzed. 





dL = -q((L«AD cosae« -L+ CD AD) sina«l 
=-yf aL + Dax | 
€ 
Б Cou + ام‎ ap» сау ЯР 
+b, 
| t 
п С, + со) h pU S eydy د کو‎ 
Су Бата ел, 52 к qu ра 
С EON = —2 (С, +С А 
Jp, 16) Cra t Co dw $e 
20 
This gives: C, -- (С +С) x 5 )5-22( 
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FIGURE 5-2 


Forces Due to Perturbation ‚u 


EIGURE SS 


Spanwise Velocity Oistribution 


Due to Yaw Rate 
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vertical tail 


FIGURE 5.4 


Vertical Tail Lever Arms 





د کا 





FIGURE 55 


Incremental Forces on Spanwise 
Element Due to the 


“Perturbation p 
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5.3 Dimensional Derivatives - Body Axes 


The analysis of the dynamic behavior of the aircraft after be- 
ing perturbed from a reference flight condition was accomplished 
through analog simulation using a body axis system. Dimensional 
parameters were used throughout the simulation. However, initial 
analysis of the variation of the non-dimensional stability derivatives 
with angle-of-attack was conducted and presented with reference to 
a Stability axis system. The method used to convert this data to 
that used for the simulation is presented in this section. 

The non-dimensional stability derivatives based on a sta- 
bility axis system were dimensionalized in a slightly different form 
than is normally used. (9) Mass and moments of inertia quantities 
were not included in order that a better physical feel for the forces 
and moments involved could be gained. The definitions of all dim- 
ensional stability derivatives used are presented in Table IIT. The 
calculation of these parmeters was programmed on a Control Data 
Corporation Digital Computer, Model 1604. The FORTRAN pro- 
gram used for this purpose is shown in Fig. 5.6 and, except for minor 
variations due to local procedures, is of standard format. (5) The 
input data consists of angle-of-attack and non-dimensional stability 
derivatives, in stability axes, inserted in the form of one-dimension- 
al, 23 element arrays corresponding to angles-of-attack from 0% to 
809 in 49 increments with angles-of-attack of 149 and 189 inserted 


in sequence. 
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The FORTRAN program of Fig. 5.6 was also designed to 
transform the dimensional derivatives based on stability axes to cor- 
responding derivatives based on body axes. The transform equations 
used were those contained in the RAeS Data Sheets and are presented 
in Table IV. (11) The angle of rotation, € , through which the axes 
were rotated was taken to be the angle-of-attack of the aircraft since 
the plane containing the MAC of the wing is within 0. 5% of coincidence 
with the plane containing the body longitudinal axis. (12) 

Tables V, VI, and VII present the dimensional stability der- 
ivatives and associated angles-of-attack in body axes for relative 
velocities of 100, 60, and 20 knots respectively. 

Several derivatives which are peculiar to the particular air- 
craft used in the investigation also require definition. These der- 
ivatives are four in number and consist of the nose fan pitch damping, 
nose and main fans lift damping, and main fans roll and yaw damping. 
These contributions were derived analytically and experimentally by 


RAC and only a brief discussion of each is presented here. (13) 


Nose Fan Pitch Damping 
From wind tunnel tests, RAC found that the nose fan pitch- 


ing moment contribution could be expressed as follows. 


S 
М. = ^C, 2 5. Ke (5-23) 


where: 
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$ 
A Cu = 6.25 ШЕ Т We Sin RAE o 


and: ( um = S 
d | Г) b/g 


Expressing о е 28 the sum of a reference flight angle-of-attack 
ToT 

and a perturbation angle-of-attack, expanding in a trigonometric 

identity, and differentiating Equation (5-23) 


with respect to e, yields: 
Ma, T6284 CUM ( - Sine, Sins, ہے‎ Cose, Cose ) (5-24) 
Using a small perturbation assumption, this may be expressed: 


ea = 6.25 4 Sue Qe ( cos SE ZA с Sin, ) (5-25) 


At low angles-of-attack, the term ie FILE IS very much smaller 
than coso Since both the perturbation, تن کے‎ , and Sin oy terms are 

small. This term may then be neglected in the low angle-of-attack 
region. At high reference angles-of-attack, the perturbation, œ 


ме? 


approaches zero. This may be seen from the expression for the 
7%. 


angle-of-attack perturbation at the nose fan due to a pitch rate on 


the aircraft as derived in Fig. 5.7 and given as: 


= а Хы COS cx, (5-26) 
М.Е. ү 


At high angles-of-attack, then, the term«x, Sind, approaches zero. 
МЕ 
Therefore, throughout the angle-of-attack range 09 « œ < 900 the 


derivative M... L may be approximated by: 
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M oce 20 % Sue. مر‎ COS en (5227) 


ше. 


The contribution of a perturbation in vertical velocity to the 
angle-of-attack change of the nose fan has been neglected in consid- 


ering pitch damping since it is very small. 


Nose and Main Fans Lift Damping 

This contribution, as well as the others to follow, was ob- 
tained by RAC directly from wind tunnel data and estimates of the 
changes of axial flow velocities through the fans. RAC shows that 


for maximum fan thrust output, the following data is obtained. 


une = 8.23 lb/Ft/sec Perfan 


4 = (d.45 Lb/ Ft [Sec 
N.F 
For the analog simulation in this analysis, a nominal total value 


is used. 


72 = — {4.0 Lb /Ft/sec 


UF eons 

Main Fans Roll Damping 
Roll damping provided by the main fans, or lift fans, results 
from the differential lift damping forces produced by the fans during 
a rolling maneuver multiplied by an appropriate moment arm. RAC 


estimated this effect from analysis of wind tunnel tests to be: 
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L Т г. 
Par. 


This value was used in the present analysis. 


Main Fans Yaw Damping 
Analysis of wind tunnel tests by RAC also yielded an estim- 


ated value of the main fans yaw damping contribution as: 
Np = - 515 Ft Lb /Rad /sec Total 


m. 


This value was likewise used in this analysis. 
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TABLE II 


Definitions of Dimensional 


Stability Derivatives 


| Longitudinal | Laterai-Directional 
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EEE 
Transform Equations 


Stability Axes to Body Axes 


XuCose - (Ху + 2 и? 9! Пе Соз е + ZwSine 
5 XwCose«(X, аше S Zy Sine 
2 созе е x Sin €C Oses XwSine 
2 Созе N + Х „) 5! пеСоз е “+ 2 
< М „Со5 е 
< 2 е Созе 
Го Созе 
2 а Cos e 


Ма 


Longitudinal 


US, 

= Ny Cose+ Ly Sine 

ES Cos E Ny SITE 

Y ues ٴ‎ = Y 5 ٣ 

= Np Cos e - (Nr - Lp)SineCose-L, Sin e 
Lp Cos 6 -(Lr * Np)Sine Cose "Nr Sine 
Yr Cos e « Yg Sine 

BN Coc € = ee 1m 
= 7 Cos e - (Np - Lp)SineCose - Np Si n'e 
C 


Lateral2 Directional 
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TABLE V 
CIMENSIQNAL STABILITY DERIVATIVES 


BCCY AXES 
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TABLE V (Continued) 
CIMENSIONAL STABILITY DERIVATIVES 
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TABLE VI 
CIMENSIONAL SIABILTIY DERIVATIVES 
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TABLE Vl (Continued) | 
CIMENSTIGNAL STABILITY DERIVATIVES 
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TABLE VI (Continued) 
CIMENSIONAL STABILITY DERIVATIVES 
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TABLE VI 
CIMENSIONAL STABILITY DERIVATIVES 
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(Continued) 
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TABLE VII (Continued) 
CIMENSICNAL STABILITY DERIVATIVES 
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5.4 Analog Simulation 


Analog simulation of the equations of motion of the aircraft 
was conducted on à Donner Analog Computer, Model 3100, in real 
time. A complete list of the equipment used is presented in Table 
VIII. Fig. 5.8 shows the computer equipment set-up for simulation 
of the longitudinal equations. 

Since the equations of motion used in the analysis are un- | 
coupled between the longitudinal and lateral-directional modes, 
simulation of the two modes was conducted separately and will be 
discussed in that manner. Techniques used in the analog circuit 


design are standard. (4) 
5.4.1 Longitudinal Mode 


For the simulation study, the longitudinal equations of 
motion, as previously presented, were rearranged and used in the 


following form. 


Ш = Хч + Xul - дсозе,е - w q (5-28) 
m m 

(m — Z5) Ur = Z u + Zur + 2 "02۰٤۷۶ (5-29) 

| = Ме» 4 МуФ + Мф + АМне (5-30) 

B X *3 в? B 


where: A Mye = ШІ ТЕ c qw. e, 


The fan contributions have been kept as separate entities 
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in order to facilitate easy removal of these quantities from the 
circuit to determine the effect of fan damping on the system. 

Figs. 5.9 and 5. 10 show the complete analog circuit for the 
longitudinal mode. For a given angle-of-attack and free stream 
velocity, the minor derivatives, Xu 2, А ad and UNT were 
assumed constant and inserted in the circuit as simple potentio- 
meter settings. The major derivatives, AZ 3 ИЕ Me, and 
Me, were not assumed constant. Their variation for + 1? of angle- 
of-attack from the reference angle-of-attack was programmed into 
five variable base, diode function generators. The derivatives 


Zt, and 2, were approximated by proper scaling of M and Ma 


respectively since, as a good approximation: 


Z: = С Mo 
eg ۲ی‎ 


E , and 2 , are constants of proportionality defined by the above 
relations, This, of course, is not an exact representation of Loe 
or La but the approximation is considerably better than assuming 
Z and Le to be constants as the aircraft deviates from the refer- 
ence flight condition. The inability to program all derivatives on 
function generators was dictated by the limited number of function 
generators available. 


The aircraft was perturbed from its reference flight 


condition by a step input or impulse of pitch attitude, © , by means 
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of the switch network shown in Fig. 5.9. 

Table IX designates the reference flight condition for each 
investigative run. Also indicated is nose and main fan damping 
status and the © perturbation used. Tables X, XI, and XII list the 
levels and potentiometer settings used in the various runs. 

In programming the function generators for the major der- 
ivatives the values of a given derivative at the reference angle-of- 
attack and + 1° from this angle were found from Tables V, VI, and 
VII of Section 5. 3 and a convenient maximum value chosen upon 
which the voltage level was calculated. This was done to permit 
the function generators to operate near, but not at, the maximum 
voltage output of 100 volts. As an example, suppose for an 89 


reference angle-of-attack a parameter, h » varied as below. 


h=|.2‏ کے بے 
ос = 8" h = I. 3‏ 
x = 42 he ls‏ 


Then the level of h , а ; would be based on a maximum value of 
h of 2.0 giving: 


Ih = loo volts = 50 volts /unıt h 
А 
The maximum voltage output from the function generator, 95 volts, 


would then occur at & = 8°, 
Since the function generators were programmed with respect 


to a specific reference angle-of-attack only the change, or pertur- 


bation, of angle-of-attack information needed to be channeled to the 
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generator as the “Y variable. This perturbation, defined in а 
manner similar to the perturbation of nose fan angle-of-attack 


shown in Section 5.3, is: 
A = ur COS cx, (5-31) 
V 
5.4.2 Lateral-Directional Mode 


The rearranged lateral-directional equations of motion 


used in the analog simulation are as follows. 


АР = (Ур eur) p + (Ye - u) r+qcose $ + q Sine, y (5-32) 


Er + Luar + P Ben ° 

PA Шығады вее 27 

г <Ер + Ммм + Мро+ Мго+ Мъри е) 
С. С С C C 


Once again the fan contributions are kept separate for the afore- 
mentioned reasons. 

Fig. 5.11 depicts the complete analog circuit for the lateral- 
directional mode. Since this mode was assumed uncoupled from 
the longitudinal mode, motion in the lateral-directional case will 
not change the aircraft angle-of-attack. For a given reference angle- 
of attack and flight velocity, then, all the lateral-directional deriva- 
tives are constant. Variation of these quantities with sideslip angle 
or roll angle is also negligible for small perturbations. Hence, all 
derivatives in the lateral-directional mode were inserted as simple 


potentiometer settings. The circuit as shown in Fig. 5.11 is 
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applicable to the areas of flight in which the derivative signs are 
identical to those at of = 8°. To accommodate changes of sign, 
the appropriate feedback lines had to be adjusted physically. 

The aircraft was perturbed from its reference flight con- 
dition by a step or impulse input of lateral velocity, V , which, 
in effect, perturbed the aircraft in sideslip angle @ : 

Table XIII lists the reference flight condition for each in- 
vestigative run. In addition, the status of fan damping and the ^7 


perturbation used is also indicated. Tables XIV and XV present 


the levels and potentiometer settings used in the various runs. 
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TABLE VIII 


LIST OF ANALOG COMPUTER EQUIPMENT 


Donner Scientific Company Analog Computer - Model 
3100 - with associated equipment: 


a. Dual DC Amplifier, model 3101 (15) 
b. Reference Supply, Model 3180 

c. Power Supply, Model 3105 

d. Power Panel, Model 3115 


Donner Variable Base Diode Function Generators - 
Models 35, 3750, and 3751 


Donner Electronic Multipliers - Models 33, 3730, 
and 3735 with associated Control Unit, Model 3736 


Data Equipment Company Plotamatic x-y plotter - 
Model 600M 


Brusch Instruments Mark 200 8-channel pen 
recorder 
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TABLE IX 


Analog Simulation 
Reference Flight Conditions 


Longitudinal Mode 


ہد کت کو جو تا 

Knots یں‎ M - A = 
j и ا اب‎ 
A 
ОСС С ЕС, E YS EOS 
ЖЕТЕН 


ei [O IN 
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—LB——— 4 ې ہت‎ 
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TABLE X 
Analog Simulation 
Invariant Levels 
Longitudinal Mode 


100 Volt/Deg 5730 Volt/ Rad 
10 Volt/Deg 573 Volt/ Rad 
ЕТ о 
Ө 100 Volt/Deg/Sec 5730 Volt/Rad/Sec 
Рд о 
100 Volt/Deg 5730 Volt/ Rad 
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TABLE Al 


Analog Simulation 


Variable Levels. 


Longitudinal Mode 
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TABLE XIL 
Potentiometer Coefficients and Settings 
and Amplifier Gains 
Analog Simulation — Longitudinal Mode 


Pot Number 


Problem Coefficient 


Pot Setting/Amplifier Gain 
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TABLE XII 


(continued) 


Pot: Number 
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TABLE XII 


(continued) 


Pot. Number 


Run Probiem Coefficient 





Pot. Seiting/Amplifier Gain 

0 08 0 2 ema. 
LH ١ ٢ ٢ 
K AI 
pos | | ول‎ ЕР 1817 || 
Na | dune ui 
EM 2-5 --1- 1-6 

fee өш вы 

и таз 7٤۲ 
Tor e ЕЕ HE Ыы 
ES ONS 
pee ei TTT sees tes 
БЕС | | ۷۳۷۳ Ji] Y ||| t [[p425 |1]|2723| 5. 
L. КЫ аре 
[pee 9 3 
tow | | || [208 |! | ||| | ||| | kore зот 
ЖЕН ЕЕЕ Pp] tT tT 
con 77 | | [1 | ۹۹ ۱ ۹ ۹ ۹۹۹۹ ۹ ۵ 
اہلاس‎ eest e 
tos ات‎ 7 TG 
не | ааа. 
ھ""ھ‎ Seer) ۲ Е ЕЯ 
Раа, 
lil ti ty Pte 


115 





ТАБЕЕ ЖЫ 


(continued) 
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TABLE ALLE 


Analog Simulation 
Reference Flight Conditions 


Lateral-Directional Mode 


Lift Fan v Pert 


Roll I Yaw Step / Imp 
: Damping Fri / Sec 


5.0 Step 
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TABLE XIV 
Invariant Levels—Analag Simulation 


Lateral Directional Mode 


7 
Lp 
" 
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TABLE XV 
Potentiometer Coefficients and Settings 


and Amplitier Gains—Analog Simulation 


Lateral-DíPectional Mode 
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TABLE XV 


(continued) 
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5,5 Analog Simulation Strip Recorder Traces 


The results of the analog simulation are shown in the follow- 
ing strip recorder traces. For the longitudinal mode, perturbations 
of Uu. ft. /sec., w- ft. /sec., o« degrees, O degrees, and q degrees/sec. 
were recorded in response to a step or impulse perturbation of ©. 
The longitudinal responses are shown in Figs. 5.12 through 5. 34 and 
are designated Run LO 1 through LO 19a. Table IX lists the reference 
flight conditions corresponding to each run. 

For the lateral-directional mode, perturbations vft./sec., 

Q degrees, P degrees/sec., Y degrees, and Гдергеев/вес. The 
lateral-directional responses are shown in Figs. 5.35 through 5,57 
and are designated Run LD 1 through LDI&. Table XIII lists the 


reference flight conditions corresponding to each run. 
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